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Abstract 

We propose a new mechanism producing a non-vanishing lepton number asym- 
metry, based on decays of heavy Majorana neutrinos. If they are produced out of 
equihbrium, as occurs in preheating scenario, and are superpositions of mass eigen- 
states rapidly decaying, their decay rates contains interference terms provided the 
mass differences Am are small compared to widths T. The resulting lepton asym- 
metry, which is the analogue of the time-integrated CP asymmetry in — 
system, is found to be proportional to Am/F. 



PACS number (s): 98.80.Cq; 98.80.-k; I4.6O.P 



1 Introduction 



A possible mechanism leading to the production, in the early universe, of a baryon- 
antibaryon asymmetry can be found in terms of lepton number production by heavy 
Majorana neutrino decays As well known, the lepton number Lq so generated is 
reprocessed by sphaleron transitions, and partially converted into baryon number B 



where Ug is the number of fermion generations, and tih denotes the number of Higgs 
doublets of electroweak Standard Model. There are two crucial issues in this mechanism 
which eventually determine its efficiency in producing the value for the baryon to photon 
number rj ~ (10^^° 10^^), fixed by observations on light nuclei abundances. 

First of all, the number density of heavy right-handed Majorana neutrinos Ni depends 
on the mechanism of reheating. Since the produced lepton number is proportional to the 
number of Ni per comoving volume, it is quite important to have a prediction for this 
parameter, as function of neutrino masses. This has been considered by many authors 
in the usual scenario of neutrino production via thermal excitations in the bath after 
reheating. In this case neutrino masses should be smaller than the maximal temperature 
obtained during the reheating of the universe. It has been recently shown by following 
in detail the reheating mechanism that a reasonable estimate for this upper bound is 
of the order of 10^ T^h, where Trh — y^T~mpi, is the so-called reheating temperature, T 
being the decay rate of the infiaton. The value of Tjih is constrained to be in the range 
(10^ -7- 10^°) GeV in order to avoid an overproduction of gravitinos in SUSY scenarios 0, 
which would imply that very heavy neutrinos, with masses larger than 10^^ GeV would 
not be significantly produced in the reheating. Even stronger constraints on T^h have 
been obtained by considering non thermal production of gravitinos [Q. 

More recently, the neutrino production has been also analyzed in the so-called pre- 
heating scenario which corresponds to a resonant particle production during the first 
infiaton oscillations around the minimum of the potential. The main result of this mech- 
anism is that even neutrinos with masses of the order of 10^^ GeV can be efficiently 
produced in highly non-thermal momentum distribution. As will be clear in the follow- 
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ing, this feature of Ni energy spectrum distribution is crucial for our purposes. 
The second issue is related to the presence of CP violating contributions to the neutrino 
decay channels in massless fermions and Higgs bosons, and to the role of interference 
effects which make this CP violation observable. There are two such contributions to the 
microscopic asymmetry e which have been considered in literature. The first one identified 
|jl[] is due to interference of the tree-level amplitudes with the absorptive part of the one 
loop vertex (e'-like effect) 

/ 2 ^ 



with mi the mass of the lightest Majorana neutrino, f{x) can be found in Ref. |l|] and h 
is the coupling of Ni to massless left-handed leptons and Higgs bosons. 

Recently, it has been observed that one should also consider contributions coming via 
interference with one loop self-energy (e-like effect) 0-0. For example for each neutrino 
A^,- one has 
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It is worth observing that both asymmetries (H) and (|^) very much resemble the CP 
violation asymmetries for charged B mesons 

^ - /) - ns- ^ /-) 

with / denoting an arbitrary final state. As for e„ and e^, in order to have a not vanishing 
result one needs more than one contribution to the exclusive decay channel /. In complete 
analogy to lepton asymmetries, the additional contribution to exclusive decay channels 
is provided by radiative processes, the so-called penguin diagrams. 

In this paper we consider a different scenario for neutrino decays, which provides al- 
ready at tree-level in the amplitudes a microscopic contribution to the lepton asymmetry. 
The mechanism, following again the analogy with B physics, is much reminding the time 
integrated CP asymmetry in the B^ — B^ system |T^. In this case, differently from the 
charged -B, by virtue of B^ — B^ oscillations one can produce CP asymmetries already at 
tree-level. 
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We consider the case in which heavy Majorana neutrinos are produced out of thermal 
equihbrium through a preheating mechanism. For arbitrary couphngs of these neutrinos 
to the inflaton field, they are produced as superpositions of the mass eigenstates Ni, that 
we denote as inflaton eigenstates, Na- If neutrino lifetimes are less than the typical de- 
coherence time due to scattering processes in the medium, these states, once produced, 
propagate as coherent superpositions of the Ni till their eventual decay, if neutrino life- 
times are less than the typical decoherence time due to scatterings in the medium. This 
constraint leads to conditions on both the couplings h and the Majorana masses which, 
however, are neither particularly severe nor fine tuned. 

As in the usual scenario, lepton number is produced via decays in (massless) left- 
handed fermions ipij {j being the family index) and Higgs bosons $, and their C-conjugate 
particles, right-handed antifermions ip'^j and giving rise to a microscopic asymmetry, 
for each a 

Since each Na is a quantum superposition of the mass eigenstates Ni, their decay ampli- 
tudes are linear combinations of A{Ni —>■ ipij^, '^rj^'^)- Provided the Yukawa matrix 
explicitly breaks CP invariance, a not vanishing value for can be obtained at tree-level 
only if these different amplitudes may interfere. 

It should be stressed at this point that it is always difficult to deal with unstable 
particles in the clean framework of Quantum Field Theory, since, in this case, they can- 
not be identified with asymptotic states in some Hilbert space, but rather they appear 
as resonances in S-matrix elements. An evaluation of the asymmetries and in this 



approach has been done in [^. In our case, one should consider the decay of the infla- 
ton fleld into the states ipLj^i '^rj^'^ which resonantly proceeds via intermediate N^- If, 
however, the heavy neutrinos have very small decay widths compared with their masses, 
which we will always assume in the following, one may envisage this process as the pro- 
duction of the (quasi)-stable state A^^ and its subsequent decay. To describe this process 
we will introduce a neutrino wave function, accounting for its exponential decay, whose 
espression can be deduced from the form of the corresponding propagator in the vicinity 
of the resonance. 



4 



The basic idea of the paper is as follows. Let us consider an arbitrary final state, like 
ipLj^ or ipj^j^'^, with invariant mass /i. It can be produced by the decay of neutrinos 
with mass mj if /i is, say, in the range rrii — Vi/2 < ^ < rrii -\- Vi/2. If there are two 
such neutrinos Ni and A^fc, whose mass difference is smaller than their average width 
^ik = (r^ + rA:)/2, it is impossible to distinguish if this decay product is the result of the 
decay of Ni or rather A^^. In other words, for the process ipLj^ one expects in this 

case the contribution of the amplitudes due to both possible processes: ^ Nj, ^ ipLj^ 
and i/^lj^. 

The interference between the two amplitudes leads, in general, to not vanishing micro- 
scopic asymmetries Ca, which are found to be proportional to the factors (mj — mj)/Tij 
and to the imaginary parts of the matrix products hh) . Of course, if the Majorana neu- 
trinos would be absolutely stable states, or, in the case of decaying particles, if their mass 
differences are much greater than their widths Fj, any interference would be impossible 
and the resulting asymmetry would vanish. 

The paper is organized as follows. In section 2 we briefly summarize the main features 
of the model describing heavy Majorana neutrino dynamics, whose out of equilibrium 
production mechanism is discussed in section 3. The estimate of the resulting macroscopic 
L asymmetry is presented in section 4. Finally, section 5 contains our conclusions and 
outlooks. 

2 The model 

Let us consider heavy neutrinos ujn and z/£j {i denotes the family index) with a Majorana 
mass term 

Lm^- {t^u Mij VRj + VRi Mij i/£^.) , (6) 

with M a symmetric real matrix. This mass term can be diagonalized in terms of a set 
of Majorana neutrinos Ni, with a latin letter as family index, defined as follows 

r^m = PRW,,N, , uI, = PlW,,Nj , (7) 

where W is an orthogonal matrix such that W'^MW is diagonal, and Pr^l = (l±75)/2. 
Denoting with x the scalar field behaving, in a certain period of the evolution of the uni- 
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verse, as the inflaton, the production of the heavy neutrinos via the reheating mechanism 
takes place due to a Yukawa term in the Lagrangian density of the form 

Lx = -X {T^li Gij ^Rj + T^m Gij J'lj) , (8) 

where again G is a symmetric real matrix. Let us denote with A^^^ (with a greek letter as 
index) the basis of G eigenstates. In this basis the Majorana mass matrix M is in general 
not diagonal, the two basis being connected by an orthogonal transformation 

N^ = U^,N, . (9) 

Finally, the heavy neutrinos are also coupled to massless left-handed leptons and to the 
standard SU (2) Higgs doublet $ through a Dirac term 

Ln = - (^L-*) K.N, - j h,,N, , (10) 

with h the Yukawa coupling matrix in the family space and 




3 Neutrino production 

In the usual scenario, neutrinos, along with all species of particles with masses below 
the maximal temperature achieved during reheating, are produced as thermal excitations 
when the inflaton releases its energy density and the radiation epoch starts. This means 
that, if the iVj have masses smaller than this temperature, they can be thermally excited. 
In this case, the Ni field configuration would correspond to a thermal distribution of 
particles with definite mass given by the eigenvalues of the Majorana mass matrix 
M. When eventually the temperature decreases down to rrii, if the decays into massless 
leptons and Higgses take place in a out-of-equilibrium condition, then a macroscopic 
lepton asymmetry can be produced, provided all Sacharov conditions are satisfied. This 
scenario has been widely studied, and, as mentioned in the previous section, the effect is 
basically due to the interference between the tree-level decay amplitudes and the one-loop 
contributions • In light of the bounds on Tr// and the discussion in the Introduction, 
this scenario seems only viable for Majorana neutrinos with masses lighter than 10^^ GeV. 
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We would rather analyze the case when these neutrinos are produced through a 
preheating mechanism as coherent superpositions of mass eigenstates. In this scenario 
particles are resonantly excited due to the oscillatory behaviour of the inflaton field 
x{t) = Xocos{m^t), being the second derivative of the inflaton potential at its mini- 
mum. In the following we will consider the simplest case of a quadratic potential for x- It 



has been shown by several authors, @], g], [[T^l^lllll! that is possible to produce spin(O), 
spin(l/2), as well as spin(3/2) particles through non-perturbative effects. The neutrino 
field satisfies a Dirac equation on a Friedmann-Robertson-Walker spacetime with an ef- 
fective time dependent mass matrix A4. In conformal time drj = dt/R{t), with R{t) the 
scale factor 



i . 3 rfi? n 



Ni^Mij(il)Ni , (12) 



with 

M,,{v) = M^J+x{v)G^, . (13) 

Equation ( ]T^ describes oscillators with a complex time varying frequency. In the simple 
case of diagonal M and G matrices the neutrino number density is obtained by using a 
time dependent Bogoliubov canonical transformation |0. In the general case, however, 
since M and G are not necessarily simultaneously diagonal one should first diagonalize 
via an orthogonal transformation the effective mass Ai, whose eigenvalues may have quite 
an involved dependence on xiv)- L^t us consider the two possible cases: 

i) if the order of magnitude of the matrix elements of M is much larger than the one 
of xiv) G-i may treat the term xiv) G as a perturbation to M . Thus in the 
Majorana mass eigenbasis, up to the first order in perturbation theory we get 

Mij{vi) c^[mi + x{ri)Gi^ 5ij . (14) 

Hence the Ug equations for i = 1, 2, .., tt.^ are decoupled and can be treated as in Ref. 
[^. In this case the inflaton would produce neutrinos already as mass eigenstates, 
but the resonant condition det [M. {r])] = for an explosive production of heavy 
neutrinos is never satisfied. 

ii) in the opposite case, when the time dependent mass term, provided by the coupling 
to the inflaton, is comparable or even dominant over M, it is convenient to use the 



inflaton basis to write the mass term as 

Ma^if]) - [MafS + Xiv) gJafs] ■ (15) 

Thus, the preheating neutrino production occurs for those values of r] for which 
det [A^(?7)] vanishes. 

In case ii), if we assume, as it will be clear in the following, small off-diagonal terms 
in the matrix M compared to the diagonal ones, we can write the preheating production 
condition as 

det[Miv)] = KUiM^a + xiv)9a) + o('^)=0 , (16) 

where m^^ just denotes the order of magnitude of the off-diagonal mass terms in M and 
m the one of the diagonal entries. At lowest order in the ratio rriais /m, eq. (|T6|) is satisfied 
if, for some r^* 

M^o. + Xiv*) 9a = . (17) 

In this case heavy Majorana neutrinos will be resonantly produced. However, since 
the conditions (|16D is only satisfied up to terms of the order map/m, this implies that the 
production rates are suppressed to some extent. In other words it is as these neutrinos 
at r/^, would not be produced as massless but with a mass of the order of m^^/m. This 
effect which leads to an exponential suppression factor of the number of heavy Majorana 
neutrinos produced 0, can be neglected if the following condition is satisfied 



« \l9a X'iV*) ■ (1^ 



m 

Since QaX'iv*) = 9cJT^xX{'n*) ~ ^ above condition becomes 



m \ m J 

The constraint (^) gives for example for the typical values = 10^^ GeV and m = 
10^^ GeV, nictis/rn < 10^^, which does not severely affects the order of magnitude of the 
off-diagonal elements in Mq,^, and still allows for a quite large mixing, given by the non 
diagonal elements of the matrix [/, see eq. (||). 

Under the condition (|T9|) one can safely apply the results of Ref. where to solve 
eq. (0), in case ii), one writes down the momentum expansion for Majorana neutrinos 
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Na{ri, x). By a time dependent Bogoliubov transformation it is then possible to diagonal- 
ize the Hamiltonian in terms of quasi-particle creation and annihilation operators and, 
with a customary procedure, to deduce the number of produced particles as the expecta- 
tion value of the particle number operator on the quasi-particle vacuum. Depending 
on the value of the parameter qa, defined as 



..-44^ , (20) 

4 m. 



X 

with x(0) the initial value of the inflaton field configuration, the final number density 
quite rapidly reaches the bound due to Pauli blocking and can be expressed in terms of 
the maximal momentum kmax of the distribution 



na ^ kl^ax - m^Mo^a ■ (21) 



At later times, when eventually oscillations are damped to smaller values, the fraction of 
the inflaton energy density transferred to heavy neutrinos, Pa/ Pxi is frozen to the value, 
see Ref. 11, 



Pa m 



2 

qa . (22) 



Px x'iO) 

In the framework of a chaotic inflation scenario, from the observed amplitude of the 
density perturbations on large scales, is constrained to be of the order of 10^^ GeV. 
Furthermore, with x(0) ~ mpi, one gets pa/ Px — 10~^^5'a. In the following this result will 
be used to estimate the final lepton number produced by heavy neutrinos. 

4 Neutrino propagation and decay 

In the previous section we have described the production mechanism of the Na neutrinos 
regarded as asymptotic stable states. In principle the neutrino production by the inflaton, 
and their subsequent decays in fermion-Higgs pairs should be considered as a whole, x ~^ 
Na '4^ Lj^ 1 '4'%^'^ ■ The metastable character of the A^^ would therefore be encoded in 
the expression of its propagator. However it should be first noted that the nonperturbative 
mechanism outlined in the previous section is basically istantaneous If in addition 
the Na are sufficiently long-lived, i.e. the widths are small compared with their masses. 
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we can safely describe the process as occurring in three stages: the production of the 
quasi-stable states Na, its propagation in the medium and its eventual decay. 

If neutrinos are produced in the inflaton basis A^^, say at time t = and with mo- 
mentum k, they start evolving as a linear superposition of the mass eigenstates Nj of 
eigenvalue rrij as follows 

Na,{x] k) = U^j Nj{x- k) , (23) 

with X = (t, x) and Nj{x] k) the wave-function of the mass eigenstate Nj. The form of this 
wavefunction, which account for the metastable behaviour of Nj and can be inferred from 
the form of the propagator, will be discussed in the following. In presence of a medium, 
scattering processes will tend to destroy the coherence among the components of the wave 
function Na{x; k) before its decay, and to populate the universe with thermalized neutrino 
mass eigenstates. In order to avoid this, one has to impose the condition 

T^»n{a)=T,, , (24) 

where Va denotes the Na decay rate 

Ta = Ua,{Y.{fm{N,\f)]jUj^ , (25) 

where the sum is over all possible final states /. With a we denote the cross section of 
the relevant scattering processes, averaged over the incoming particles distribution with 
number density n. The dominant contribution to Fq, comes from the two body decay 
channels ipu^- If the neutrino masses are of the same order of magnitude irii c:^fn'^i, the 
factor in bracket in eq. (^) simplifies to 

J2{fm{N,\f) ^ {hh% - ^ r„ ^ {hh^u ^ ■ (26) 

Notice that the case of almost mass degenerate neutrinos is actually the scenario we are 
mostly interested in this paper. 

The main contributions to Tsc correspond to the processes N^ $ — > Ni $, N^ipLi — ^ 
Ni ipLj and crossed channels, as well as to scattering over the large amount of inflaton 



quanta x- Condition (|25D constraints more severely the Yukawa couplings h and G when 



the mean energy of the massless ^, ipii and x is larger than the neutrino masses. This is 
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due to the rapid increasing with this mean energy of the product ^#,V)l,x ('^)- isuct, defin- 
ing, in the Na rest frame, the effective temperature T*, which represents the mean energy 
of the massless fermions or Higgses, and similarly the one of the infiaton excitations, 
we have, for scattering over massless fermions and Higgses 

1 



Stt m 

and similarly, for scattering processes over x bosons 



(27) 



. (28) 

For low momentum neutrinos and also give the mean energy in the comoving frame. 
Actually from the discussion of section 3, we see that in our preferred choice for m > 10^^ 
GeV and ~ 10^^ GeV, neutrinos are basically emitted as non-relativistic particles, 
since kmax/^ — {m^jmY'^'^ < 1, thus we can safely neglect any effect due to the difference 
between the neutrino rest frame and the comoving frame. Since the number density of 
incoming particles can be expressed as n^^ ~ r;,$ ~ T^^ ~ T^, the condition for no 
decoherence becomes 

No particularly fine tuned condition on the Yukawa couplings h and G follows from 
(p9|) even for very heavy Majorana neutrinos, m ~ 10^^ GeV. Assuming in fact the 
extreme limit ~ (10^^ ^ 10^^) mpi which, still compatible with a classical description 
of spacetime structure, seems to be suggested by an exponential production of x quanta by 
the infiaton oscillating configuration |T^, we get, as an extremely conservative estimate, 
that there is no decoherence if we take the Yukawa couplings G, h up to (10^'^ ^ 10^^). 
This implies the conservative upper bound for (0) 

< (lO^ - 10«) . (30) 

If condition (p9D is satisfied, heavy neutrino states (^3]) evolve coherently till their decays 
into massless fermions and Higgses. 

To evaluate the decay rates in these channels we first have to specify the expression for 
the wave functions Nj{x] k). For this purpose we first recall that, in a neighborhood of the 
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mass eigenvalues nij, the Nj propagator Gjij?) gets the familiar Breit-Wigner behaviour. 
Using Lehman spectral decomposition we can also write Gj{p'^) as the superposition of 
free propagators of given masses via two spectral functions 

gap') - , ^tT- r = r ^1^^') + ^?(^')) 2 , (31) 

— rrij + imjLj Jo \ J y _ _l_ 

with C the charge conjugation operator. For /x^ ~ m^, it is straightforward to get 



r ■ (32) 



This result follows by expanding the free propagators as the principal value of {p^ — IJ^)~^ 
and its Dirac 5-contribution. 

The positive definite spectral function can be simply related to the wave function 
Pj of the unstable state Nj in configuration space. If we consider the expansion of the 
field Mj, which represents the Major ana neutrinos, as follows 

•^^■^"^^ I^^'^"^' + ^^^^ 

with Nj{x;k) a set of neutrino wavefunction with spatial momentum k, encoding its 
unstable character, it is straightforward to see that we recover the correct Breit-Wigner 
behaviour of the corresponding propagator if we represent these wavefunctions as 

Njix;k)^ Jdp'pj{p')ule-''-^ , (34) 

with |pj(//^)p = Wj{p^). We have denoted with u'i{x) a four component spinor solution 

— * 

of the Dirac equation in momentum space with mass p and momentum k. The spectral 
function wj does not uniquely fix the wavefunctions Nj{x;k). The simplest choice for 
Pj{p^), satisfying Ipjip"^)]^ = Wj{p'^), is again a Breit-Wigner function 

- ^^^^H:^-TT • (35) 

We stress that in what follows it is not really crucial the particular choice for pjlp"^), 
satisfying = Wj. The only relevant aspect we have to require in order to get a non 
vanishing lepton asymmetry is that Pj{p^) has an imaginary part depending on the two 
parameters ruj and Tj. This assumption is quite robust, if we one recall the fact that phase 
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shifts due to the production of any resonance is strongly dependent on these quantities. 
Finally we observe that in the limit of vanishing decay widths Pj(/i^) —>■ — so in 
this limit Nj{x] k) — > m^-' (x)e~*'^'^, which describes a stable neutrino with mass rrij and 
momentum k 

We can now evaluate the decay rates of the N^. As we have discussed in the preheating 
mechanism neutrinos are basically produced with low momenta, much smaller than their 
masses. It is therefore a good approximation to consider the case of neutrinos emitted at 
rest in the comoving frame. From eq. (^TJ) all corrections to the results below are at most 
of the order of (m/m^)^/^. 

Depending on the Majorana neutrino mass spectrum, it is possible that two or more 
Pj significantly overlap. This occurs whenever \mi — mj\ < (Fj + Fj)/2. Using (0), the 
total decay rate of Na into pairs ipLp^ as well as into the C-conjugated channels 
is given by 

rig Ug 

J2TiN^^ij,,^)=Ei^h, , J2nNc. ^ ^Rp^'^) = {Ei^y , (36) 

p=l p=l 

where 

2q = Uai{hh^)ijUaj = {^S] ) o"^^^ ^ and j , (37) 



In eq. ( |5^ ) the factor /i /(Stta/ |A;P + jj?) ~ /^/(Svr) is the result of integration over the 



decreasing functions for fi"^ << mf^rn?^ and ji^ » ,,o^, 



phase space for final massless particles with fixed initial mass fi and momentum k. Notice 
also that we have integrated in lij over the range [0, cxd[, rather than in the narrow 
neighborhood of the masses rrii and rrij. This is justified by observing that 1^ receive the 
main contribution from these neighborhoods, since the products pi{p^)p*{p^) are rapidly 

2 2 

A straightforward computation shows that the microscopic asymmetries are given 

by 

^ 2E.<,Xm[S'i]Xm[J,,] 

2 E.<, Tie [J,,] + Y.Zi '^l hi 

From this result we get that, in order to have at tree-level not vanishing microscopic 

asymmetries it is necessary that: 
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i) a. CP violating term is contained in the Yukawa couplings h; 

ii) the integrals lij of the spectral functions contain a not vanishing imaginary part. 

As already stated, this second condition is realized if at least two of the neutrino masses 
satisfy the condition \mj — mi\ < (rj + rj)/2, otherwise the two kernels Piifi"^) PjifJ''^) have 
no significant overlap. If rrij > rrii, the main contribution to the imaginary part of 1^ is 
expected for rrii < p < rrij. In this interval the phase difference of pi and pj is almost vr. 
Furthermore, even if rrii = rrij, an imaginary part for lij is expected if the two widths are 
sensibly different. 

Defining = (Fj + Tj)/2, rriij = {rrii + rnj)/2, A,j = {mj - mi)/(Ti + Tj) and 
jij = (Tj — Ti)/Tij, and if we also assume for simplicity that both Ajj, •jij < 1, the 
asymmetry can be obtained as an expansion in powers of these parameters. In the narrow 
width limit, a simple calculation up to the first order in Aij and •jij gives 

It is interesting to consider the case of only two generations, for which the expression of 
the L microscopic asymmetry is particularly simple and the orthogonal matrix U is given 
by 

TT _ ( cos{9) sm{9) \ , 

[- sm{e) cos{e) J • ^ ^ 

A simple calculation gives at the lowest order in A12 and 712 



2sin(2^)Jm \ihh^)u 

ei,2 = 



A12 + 7i2 



(/l/lt)n + (/^/^t)22Tsin(2e)7^e[(/^/^t)l2] V '"m^ 

± sm 2^ Acp A12 + 712 ^ ^ . . , 42 

V mi2/ [hh^)ii + [hhT )22 =F sm(26')7^e [{hh^ )i2\ 



with 



2Im 

^CP 
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{hh^)u + {hh^)22 ' ^^^^ 

representing the strength of the CP violating effects. Notice that since we are working 
with Majorana neutrinos, even for two generations a phase in the matrix h cannot be 
washed away by a simple redefinition of fields, so {hh'^)i2 is in general a complex quantity. 

There are several features of eq. (^) which is worth observing. First of all the two 
asymmetries tend to cancel each other, simply because the two numerators of ei and €2 
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are opposite by the orthogonality of the mixing matrix U. This is exactly analogous to 



GIM mechanism |rB|. However, already at first order the sum ei + 62 does not vanish 
because the denominator in eq. ( ^21) is different for the two neutrinos. Furthermore, 
the total lepton asymmetry L is given by the microscopic asymmetries weighted by the 
corresponding neutrino number densities n^. As long as rii 7^ n2 the value of L is not 
expected to vanish, though a partial cancellation still takes place. It should also be pointed 
out that at lower order there is also a contribution from the widths difference parameter 
712, which is however suppressed by the small factor Tu/mu- 

Finally notice that the value for the asymmetry is obtained at tree-level in decay 
amplitudes, while both contributions previously considered, and e^, depend on matrix 
elements of [hh'^Y in the numerator, see eq.s This means that the asymmetries 



(^) can be in principle quite large since they are not suppressed by higher powers of 
the Yukawa couplings, though there is a certain cancellation among them. Also notice 
that, differently than eq. (^, ei^2 are only linearly dependent on A12, so they are less 
suppressed in the limit of small A12. Incidentally, in the usual scheme the Yukawa matrix 
elements are constrained to be quite small (/i/i^)m < mi/mpi to have an out-of-equilibrium 
decay. This bound is no more necessary if the neutrino are already produced, as in the 
preheating scenario, in a non thermal way, and the reheating temperature Tptn is low 
enough to prevent from a subsequent thermal production of Ni. 

Using the results of section 3, we can finally estimate the total lepton number n^, 
normalized to specific entropy. Since the energy fractions pa remain constant till the 
inflaton decay into radiation, after the reheating stage one gets 

— - = — 2^ Pa ^ 10 = 2^ ea ga ■ 44 

■5 Px Q^i VlO^'^Gel// \ m ) 

Choosing in the interval ~ (10^ ^ 10^), a reheating temperature ~ 10^°Gel^ 
and for heavy right-handed Majorana neutrinos, m ~ 10^^ GeV , one gets 

rig 

E^"~(lO-^-l) ' (45) 

a=l 

where the value for the ratio ns/s ^ ul/s ~ (10^^^ ^ 10^^°) given by primordial nucle- 
osynthesis has been used. Note that the lower bound for YlaLi of eq. (^51) strongly 
decreases if one slightly releases the very conservative condition for no decoherence (|30D . 
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Moreover, since very heavy Majorana neutrinos with mass up to 10^^ GeV are still com- 
patible with preheating scenario described in Ref. it is interesting to observe that just 
increasing the value of m of one order of magnitude, by virtue of (0), (^) and (^i]) one 
reduces the lower bound for J2a=i of the same amount. 

From the above considerations it follows that the required value for the ratio ub/s 
can be easily obtained for a wide range of the involved parameters, without imposing fine 
tuned conditions. In particular it is worth noticing that no particular mass degeneracy is 
necessary, and CP violating effects of the same order of those predicted in the electroweak 
Standard Model, and measured in K'^ — system, are already sufficent to produce a 
baryon asymmetry of the correct order of magnitude. 

5 Conclusions 

In this paper we have considered a new scenario for the production of a primordial lepton 
number, based on decays of oscillating heavy Majorana neutrinos. In the framework of the 
preheating mechanism for a non-thermal production of massive fermions, we have stressed 
the possible role in leptogenesis of coherent superpositions of unstable mass eigenstate 
neutrinos. The mechanism is similar to the way an observable CP asymmetry is produced 
in the neutral B meson system, due to the 5° — 5^ oscillation in time. In fact if the 
decaying neutrinos are linear superpositions of mass eigenstates Ni, the observability of 
the CP violation is achieved if the different A^^ may interfere. 

We have shown that if at least two neutrinos have small mass difference Am, compared 
with their corresponding decay widths F, a microscopic asymmetry can be obtained at 
tree level in the Yukawa matrices, coupling these states to massless fermions and Higgses. 
The microscopic asymmetry is found to be proportional, for small mass differences to the 
ratio Am/F. On the other hand, if the neutrino states have masses quite well separated, 
with respect to their decay widths, the interference effects we have described is vanishing 
and to get a non-vanishing microscopic asymmetry one should consider interference of 
the decay amplitudes at tree level with higher loop contributions 

In this scenario, a crucial feature is that neutrinos are produced with a non-thermal 
distribution, thus the Sacharov out-of-equilibrium condition is implemented from the very 
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beginning. If, as we have considered, neutrinos rapidly decay before any thermalization 
may occur, it is possible to avoid any wash out of the final lepton number so produced 
by inverse processes or scattering, provided the maximal temperature achieved during 
reheating is smaller than the mass of the lightest of the heavy Majorana. 

It is finally worth stressing that, even using very conservative bounds on the couplings 
involved and Majorana masses, the order of magnitude of the microscopic asymmetries 
results compatible with CP violation effects in the Standard Model. Actually this result 
holds for choices of Am which are not particularly fine tuned. 
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Note added in proof 

Shortly after completing the paper, the role of neutrino instability in mixing phenomenon 
have been presented in a recent manuscript . 
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